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ABSTRACT 


Due  to  nonlinearity  in  the  eoupling,  one  of  the  vibrational  modes  of  a  straight 
symmetrie  triatomie  moleeule  ean  be  unstable  for  amplitudes  greater  than  a  threshold 
value.  The  instability  is  due  to  the  mode  parametrieally  exeiting  another  mode.  The 
threshold  amplitude  deereases  if  the  differenee  of  the  frequency  of  the  two  modes  is 
reduced.  We  consider  the  simplest  case  of  a  symmetric  rectilinear  molecule  where  the 
coupling  has  a  cubic  nonlinearity  in  addition  to  a  linear  restoring  force.  Approximate 
analytical  results  are  in  good  agreement  with  numerical  simulations  of  the  exact 
equations  of  motion,  although  in  some  cases  the  actual  behavior  fundamentally  deviates 
from  the  perturbative  theory.  Two  physical  demonstrations  of  the  instability  are 
described,  where  the  apparatus  are  a  system  of  gliders  coupled  by  springs  and  magnets  on 
an  air  track.  Possible  quantum  mechanical  implications  are  discussed.  This  work  is  a 
fundamental  generalization  of  the  parametric  instability  of  two  linearly  coupled  nonlinear 
oscillators  that  was  reported  in  a  previous  investigation. 
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I.  INTRODUCTION 


The  straight  (often  ealled  linear)  symmetrie  triatomie  moleeule  is  a  well  known 
and  well  understood  system  that  is  used  as  an  example  in  intermediate  (Fowles)(Marion) 
and  graduate-level  (Goldstein)(Landau)  meehanies  textbooks  in  determination  of  the 
normal  longitudinal  modes  of  a  system,  see  Fig.  1.  This  system  is  known  to  have  three 
longitudinal  modes  that  we  refer  to  as  the  zero-frequeney  mode,  lower-frequeney  mode 
and  upper-frequeney  mode.  The  zero-frequeney  mode  eorresponds  to  the  rigid 
translation  of  the  moleeule  along  the  longitudinal  axis  with  eonstant  veloeity.  The  lower- 
frequeney  mode  eorresponds  to  the  eenter  mass  remaining  stationary  while  the  end 
masses  oseillate  in  opposite  direetions  with  the  same  amplitude  and  angular  frequeney  (Oo 
=  (k/m)  .  The  upper-frequeney  mode  eorresponds  to  the  outer  masses  oseillating  in 
phase  and  with  equal  amplitudes  while  the  eenter  mass  oseillates  in  antiphase  with 
angular  frequeney  coi  =  (1  +  2m7M)  (Oo  and  amplitude  sealed  by  the  faetor  2m7M.  This 
eorresponds  to  zero  displaeement  of  the  eenter  mass. 

Previous  work  dealt  with  a  system  of  two  linearly  eoupled  nonlinear  oseillators, 
for  whieh  it  was  shown  that  one  of  the  two  normal  modes  is  subjeet  to  the  parametrie 
instability.  (Denardo)  It  was  thought  that  a  mode  of  a  system  would  be  subjeet  to  the 
prineipal  parametrie  instability  only  if  its  frequeney  were  twiee  the  frequeney  of  another 
mode.  The  investigation  opened  the  possibility  for  parametrie  instabilities  under 
eonditions  that  may  more  readily  oeeur  beeause  the  instability  is  favored  (i.e.,  the 
amplitude  threshold  is  lower)  when  the  frequeneies  of  the  two  modes  are  near  eaeh  other. 
This  instability  may  have  eonsequenees  regarding  aetual  moleeules  (see  below). 

In  eontrast  to  previous  work,  this  system  does  not  deal  with  eoupled  oseillators. 
Rather,  the  nonlinearity  exists  in  the  eoupling  of  the  masses.  The  overall  eommon 
eharaeteristie  between  this  system  and  the  previous  one  is  that  the  nonlinearity  eauses  the 
normal  modes  to  be  eoupled.  Additionally,  the  instability  is  favored  when  the 
frequeneies  of  the  two  normal  modes  are  near  eaeh  other,  whieh  oeeurs  for  M  »  m  in  the 
triatomie  moleeule. 
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Figure  1 .  Symmetric  rectilinear  triatomic  molecule,  with  outer  masses  m,  central  mass  M,  and  spring 

constants  k:  (a.)  translational  motion  (b.)  lower-frequency  motion  (c.)  upper-frequency  motion. 

This  thesis  shows  that  the  nonlinearity  in  the  eoupling  eauses  one  of  the  two 
oseillatory  normal  modes  to  be  unstable  for  amplitudes  above  a  threshold  value.  The 
instability  leads  to  a  beating  motion  where  energy  is  transferred  baek  and  forth  between 
the  lower  and  upper  frequeney  modes.  For  simplieity,  only  a  eubie  nonlinearity  in  the 
displaeement  of  eaeh  spring  from  equilibrium  is  assumed.  The  instability  eorresponds  to 
the  exeited  mode  parametrieally  driving  the  other  mode,  whieh  is  referred  to  in  general  as 
autoparametric  resonance.  The  instability  is  investigated  analytieally  and  numerieally 
and  two  demonstration  apparatus  that  exhibit  the  effeet  are  deseribed. 

The  motivation  for  this  work  is  two-fold.  First,  the  instability  is  an  example  of 
how  nonlinearity  in  a  system  ean  lead  to  fundamentally  different  behavior.  The 
instability  is  espeeially  important  in  the  ease  of  a  elassieal  triatomie  moleeule,  beeause 
this  is  a  model  system  that  is  studied  by  nearly  all  physies  students.  Seeond,  the 
existenee  of  the  (elassieal)  instability  raises  interesting  questions  regarding  the  quantum 

meehanieal  behavior  of  an  aetual  triatomie  moleeule,  or,  in  general,  any  moleeule  that  has 

2 


more  than  one  vibrational  mode.  It  is  well  known  that  the  vibrational  motion  of 
moleeules  is  nonlinear  (Beiser),  so  the  normal  modes  will  be  eoupled.  Although  there  is 
mueh  literature  on  a  broad  variety  of  parametrie  exeitation  in  quantum  meehanieal 
systems  by  external  sourees,  (Quantum)  we  have  not  found  any  work  on  autoparametrie 
resonanee. 

As  a  speeulative  example  of  a  possible  eonsequenee  for  an  aetual  triatomie 
moleeule,  eonsider  a  gas  of  sueh  moleeules.  The  unstable  mode  may  not  fully  partieipate 
in  the  partition  of  energy  among  the  degrees  of  freedom  of  the  gas,  so  the  speeifie  heat 
would  be  less  than  the  predieted  value.  This  will  oeeur  if  the  temperature  is  suffieiently 
large  sueh  that  average  energy  of  the  mode  is  roughly  equal  to  or  greater  than  the 
threshold  value.  For  a  eubie  nonlinearity,  we  ealeulate  the  elassieal  amplitude  threshold 
in  Sees.  II.B  and  II. C. 

As  another  speeulative  example,  eonsider  the  ease  in  whieh  the  ground  state  is 
elassieally  unstable.  As  shown  in  See.  II.A,  this  is  possible  for  a  softening  eubie 
nonlinearity  if  M  is  suffieiently  greater  than  m  sueh  that  the  amplitude  threshold  is  less 
than  the  amplitude  eorresponding  to  the  ground  state  energy.  The  instability  eannot 
oeeur,  beeause  energy  eonservation  would  be  violated.  However,  the  elassieal  instability 
may  still  have  some  manifestation,  beeause  it  must  oeeur  at  higher  energies  as  a  result  of 
the  eorrespondenee  prineiple.  An  example  of  a  possible  eonsequenee  is  that,  in  the 
presenee  of  photons  or  eollisions,  there  may  be  an  enhaneement  of  the  transition 
probability  from  the  ground  state  to  a  state  of  the  other  mode. 

We  were  initially  motivated  to  investigate  the  parametrie  instability  of  a 
vibrational  mode  of  a  moleeule  beeause  we  thought  that  this  eould  lead  to  a  new  type  of 
laser;  speeifieally,  where  a  population  inversion  results  fundamentally  from  the 
instability.  However,  the  transition  probability  for  absorption  from  one  state  to  an 
exeited  state  will  still  equal  the  transition  probability  for  stimulated  emission  from  the 
exeited  state  to  the  original  state,  so  it  does  not  appear  that  the  instability  itself  ean  give 
rise  to  a  population  inversion.  However,  an  enhaneed  transition  probability  due  to  the 
instability  eould  be  benefieial  in  some  lasers  beeause  it  will  eause  the  population 
inversion  to  oeeur  more  rapidly. 
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II.  THEORY 


A,  MODEL  SYSTEM  AND  EQUATIONS  OF  MOTION 

We  consider  longitudinal  motion  of  the  system  shown  in  Fig.  1.  The  force,  F,  of 
each  spring  is  generalized  to  include  a  cubic  nonlinearity: 

F  =  -kx  +  amx^  (1) 

where  x  is  the  displacement  of  the  spring  from  equilibrium,  k  is  the  spring  constant,  m  is 
the  mass  of  the  outer  masses  and  where  the  nonlinear  coefficient  a  can  be  positive 
(“softening”),  negative  (“hardening”),  or  zero  (linear).  The  factor  of  m  in  the  cubic  term 
is  inserted  only  for  mathematical  convenience.  The  coupling  force  in  actual  molecules 
can  be  approximated  by  a  sum  of  quadratic  and  cubic  nonlinearities.  We  dealt  with  the 
simplest  case,  cubic  nonlinearity. 

The  displacements  of  the  masses  are  measured  from  points  corresponding  to 
equilibrium  of  the  system.  The  equations  of  motion  for  the  displacements  xi  and  X2  of 
the  outer  masses  are 


Xi  =-(0o(Xi-X3)  +  a(Xi-X3f  (2) 

X2  =-®o(x2-X3)  +  «(x2-X3f  (3) 

where  the  overdots  denote  differentiation  with  respect  to  time,  and  where  (Oo  =  k/m. 
Because  a  possible  uniform  motion  of  the  center  of  mass  (zero-frequency  normal  mode) 
is  irrelevant  here,  we  employ  the  center-of-mass  frame  of  reference.  This  reference 
frame  reduces  the  number  of  degrees  of  freedom  from  three  to  two  and  the  relationship 
among  the  displacements  is  m(xi  +  X2)  +  MX3  =  0,  so  X3  =  -p(xi  +  X2),  where  p  =  m/M. 
Substituting  this  expression  for  X3  into  Eqs.  (2)  and  (3)  yields 
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(4) 


= -a)^[(1  +  p)x^ +px2  +a[(1  +  p)x^ +PX2  ^ 


X2=-a)l  pXi  +  (1  +  p)x2  +  a[pXi  +  (1  +  p)x2  ^  (5) 


Transforming  Eqs.  (4)  and  (5)  to  the  linear  normal  mode  eoordinates  ^  =  (xi  -  X2)/2 
(lower-frequeney  mode)  and  r)  =  (xi  +  X2)/2  (upper-frequeney  mode)  yields,  after  some 
simple  algebra, 


'i  +  m%  =  a(i’  +  3y%yf)  (6) 

Ti  +  0)fll  =  (7) 


2  2 

where  (Oi  =  yWo  ,  and  where  y  is  a  mathematieally  eonvenient  parameter  that  speeifies 
the  ratio  of  the  masses:  Y=l+2p=l  +  2m/M,  which  ranges  from  y  =  1  for  M  =  0°  to  y  = 
00  for  M  =  0.  The  first  limit  corresponds  to  uncoupled  oscillations  of  the  m  masses,  while 
the  second  corresponds  to  the  absence  of  M  and  thus  a  reduction  to  a  single  degree  of 
freedom  of  vibration. 

The  square  frequency  of  small-amplitude  oscillations  of  the  lower-frequency 
mode  (normal  coordinate  ^)  is  (Oo^  =  k/m,  while  that  of  the  upper-frequency  mode 
(normal  coordinate  rj)  is  (Oi  =  ycOo  •  The  equations  of  motion,  Eqs.  (6)  and  (7),  for  the 
normal  modes  are  similar  although  not  identical  to  those  for  the  system  of  two  nonlinear 
oscillators  that  are  linearly  coupled.  (Denardo)  The  coupling  stiffness  between  the  two 
masses  there  plays  the  role  of  mass  ratio  p  =  m/M  =  (y  -  l)/2  here. 
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B,  LOWER-FREQUENCY  NORMAL  MODE 

The  analysis  of  the  stability  of  each  normal  mode  is  similar  to  that  in  Ref. 
(Denardo).  However,  because  the  current  case  turns  out  to  be  more  complicated,  it  is 
useful  to  alter  the  approach  by  first  performing  the  graphical  analysis  and  then  the 
analytical  analysis. 

For  the  lower- frequency  (^)  mode,  we  consider  weakly  nonlinear  oscillations  ^  = 
A  cos(a)t),  where  O)  and  A  are  constants.  The  relationship  between  the  amplitude  and 
frequency  of  these  oscillations  is  approximately  (Denardo) 


0) 


2 


(8) 


Substituting  ^  =  A  cos(a)t)  into  the  r)  equation  of  motion,  Eq.  (7),  using  the 

2  3 

identity  cos  (cot)  =  [1  +  cos(2a)t)]/2,  and  neglecting  the  r)  term  because  the  amplitude  of 
the  ri  motion  is  assumed  to  be  small,  gives 


r\  + 


3  3 

— yocA^  -  —  yocA^  cos  ( 2a)t ) 


ri  =  0. 


(9) 


This  has  the  form  of  the  Mathieu  equation  fj  +  +  acos(2a)t)Jri  =  0 ,  for 

which  it  is  known  that  parametric  excitation  occurs  if  the  approximate  threshold 
condition  |a|  >  2|f2  -  O)  |  is  met.  (Denardo)  A  convenient  graphical  representation  of  this 

is  to  consider  the  drive  plane  of  a  vs.  0)  .  Parametric  excitation  occurs  for  all  points 

2  2 

inside  the  V  formed  by  a  =  ±2(f2  -  O)  )  in  the  upper-half  plane.  The  vertex  of  the  V  has 

2 

ordinate-intercept  and  the  lines  have  slope  ±2,  figure  2. 
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Figure  2.  Graphical  representation  of  the  Mathieu  equation  for  which  the  area  inside  the  V  represents  the 

region  for  which  parametric  excitation  occurs. 


2  2 

The  parametric  threshold  condition  |a|  >  2|t2  -  O)  |  corresponding  to  Eq.  (9)  is 


3 

4 


Y  a  > 


(Ot 


3  A  2  2 

— yoA  -0) 


(10) 


As  in  the  case  in  Ref.  (Denardo),  an  interesting  fact  is  that  Eq.  (9)  reveals  that  the 
frequency  is  not  independent  of  the  drive  amplitude  a.  That  is,  the  linear  frequency  of 
the  ri  mode  is  altered  by  the  finite  amplitude  oscillations  of  the  ^  mode.  This  has  the 
effect  of  uniformly  shearing  the  V  in  the  plane  of  |a|  vs.  0)  . 

As  in  Ref.  (Denardo),  we  conveniently  unify  the  cases  of  softening  (a  >  0)  and 

2  2 

hardening  (a  <  0)  by  considering  the  plane  of  aA  vs.  to  ,  where  the  ordinate  can  be 
positive  or  negative.  The  vertex  of  the  sheared  double  V  in  the  plane  has  ordinate- 
intercept  (Oi  .  The  slopes  of  two  lines  of  the  double  V  can  be  determined  from  Eq.  (10) 
by  replacing  the  inequality  with  an  equality,  removing  the  two  absolute  values,  and 
applying  a  ±  sign  to  either  side  of  the  equation.  The  resultant  values  of  the  slopes  are 

^/9y  and  -T/Sy.  A  qualitative  graph  of  the  sheared  double  V  is  shown  in  Eig.  3. 
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aA^ 


Figure  3.  Qualitative  sketch  (not  to  scale)  of  the  amplitude-frequency  plane  for  parametric  excitation  of  the 
upper-frequency  (p)  mode  by  the  lower-frequency  (^)  mode  with  amplitude  A,  where  y  =  1  +  2m/M.  The 
excitation,  and  thus  instability  of  the  ^  mode,  occurs  for  points  on  the  heavy  line  that  lie  in  the  hatched 

region. 


The  ^  mode  aets  as  a  parametrie  drive  on  the  r)  mode.  The  possible  values  of  aA^ 
and  0)  are  not  independent,  but  are  related  by  the  amplitude-dependent  frequeney,  Eq. 
(8),  of  the  ^  mode.  This  eorresponds  to  the  line  with  ordinate-intereept  and  slope 

-All)  (Fig.  3).  The  slope  is  steeper  than  the  value  -4/3y  of  the  steeper  of  the  two  V  lines, 
beeause  y  >  1-  A  hardening  (a  <  0)  system  is  not  subjeet  to  the  instability,  beeause  the 
amplitude-dependent  frequeney  line  does  not  interseet  the  hatehed  region  of  the  V.  The 
softening  results  are  in  eontrast  to  the  previous  ease  of  two  linearly  eoupled  nonlinear 
oseillators,  (Denardo)  where  the  amplitude-dependent  frequeney  and  the  steeper  of  the 
two  lines  of  the  V  have  the  same  slope,  so  there  is  a  single  threshold.  In  the  eurrent  ease, 
there  is  potentially  both  an  upper  and  lower  threshold.  The  instability  is  predieted  to 
oeeur  only  for  amplitudes  between  the  two  threshold  values. 

To  ealeulate  the  two  thresholds  of  the  instability,  we  substitute  the  amplitude- 
dependent  frequeney,  Eq.  (8),  into  the  parametrie  exeitation  eondition,  Eq.  (10),  replaee 
the  inequality  with  an  equality,  remove  the  two  absolute  values,  and  apply  a  ±  sign  to 
either  side  of  the  equation.  The  upper  threshold  value  of  the  ordinate  in  Fig.  3  is  then 
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2  2 

found  to  be  aA  =  4a)o  /3.  According  to  Eq.  (8),  this  corresponds  to  zero  frequency,  as 
shown  in  Fig.  3.  Hence,  according  to  our  theory,  there  is  a  single  threshold  represented 
by  the  lower  value.  Parametric  excitation  of  the  r)  mode  occurs  in  the  softening  (a  >  0) 
case  when  the  amplitude  of  the  ^  mode  corresponds  to  a  value  above  the  dashed  line  in 
Fig.  3.  The  calculation  of  the  lower  threshold  yields 


aA^ 


> 


4 

3 


Y-1 

(3y-1) 


a)„ 


4  P 
3(1  +  3p) 


a)„ 


(11) 


where  we  have  recast  the  result  in  terms  of  the  ratio  of  the  masses  p  =  m/M  =  (y  -  l)/2, 
which  is  more  physically  accessible.  Because  the  energy  of  the  ^  mode  must  decrease  if 
the  ri  mode  is  excited,  Eq.  (11)  represents  the  amplitude  threshold  for  instability  of  the  ^ 
mode.  The  theory  predicts  that  the  instability  only  occurs  for  a  softening  (a  >  0)  system. 

It  should  be  noted  that  our  calculations  are  based  only  on  a  leading-order 
perturbation  theory.  The  lines  in  Fig.  3  are  actually  approximations  to  curves,  so  an 
upper  threshold  may  exist  for  O)  >  0.  The  existence  can  be  resolved  by  a  higher-order 
perturbation  calculation.  An  alternative,  which  we  pursue  in  Ch.  Ill,  is  to  perform 
numerical  simulations  of  the  exact  equations  of  motion,  Eqs.  (6)  and  (7). 

C.  UPPER-FREQUENCY  NORMAL  MODE 

The  analysis  of  the  stability  of  the  upper-frequency  (rj)  mode  proceeds  similarly 
to  the  ^  mode  in  the  previous  section.  However,  in  contrast  to  the  system  in  Ref. 
(Denardo),  the  amplitude-dependent  frequency  of  this  mode  is  not  the  same  as  that  for  the 
^  mode.  Rather,  for  weakly  nonlinear  oscillations  r)  =  Acos(a)t),  the  equation  of  motion, 
Eq.  (7),  yields 


(12) 
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which  differs  from  that  of  the  ^  mode  by  a  factor  of  in  the  amplitude  term,  so  the 
frequency  shift  is  greater  beeause  y  >  1-  This  is  due  to  the  fact  that,  for  a  specified 
amplitude  A  of  the  two  outer  masses  in  the  system  (Fig.  1),  the  springs  are  stretched  and 
compressed  a  greater  amount  for  the  upper-frequency  mode  than  the  lower-frequeney 
mode,  beeause  M  oscillates  in  antiphase  in  the  first  case  and  is  motionless  in  the  seeond 
ease.  It  is  easily  shown  that  the  amplitude  of  the  spring  oscillations  is  simply  yA.  The 
slope  of  the  line  of  the  amplitude-dependent  frequency,  Eq.  (12),  in  the  plane  of  aA  vs. 
0)^  is  -4/3y^ 

The  equation  analogous  to  Eq.  (9)  is 


a)„ 


^-y^ocA^  — y^ocA^  cos  (2a)t) 


^  =  0 


(13) 


This  is  identieal  to  Eq.  (9)  if  we  let  r)  ^  (Oi  ^  (Oo,  and  y  ^  y^  in  Eq.  (9).  Henee,  the 
slopes  of  the  sheared  double  V,  whieh  were  -4/3y  and  -4/9y  in  the  previous  ease,  must 
now  be  -4/3y^  and  -4/9y^.  Eurthermore,  by  the  same  replaeements,  the  equation 
analogous  to  Eq.  (10)  must  be 


^y^|a|A^  > 


2  3  2  A  2  2 

(Oq  -  ocA  -  0) 


(14) 


The  slope  -4/3y  of  the  amplitude-dependent  frequency  line  of  the  r)  mode  from 

Eq.  (12)  is  less  steep  than  the  slope  -4/3y  of  the  steeper  line  of  the  V,  beeause  y  >  1. 

2 

Eurthermore,  the  slope  of  the  amplitude-dependent  frequeney  line  equals  the  slope  -4/9y 
of  the  less  steep  line  of  the  V  when  y  =  3.  There  are  thus  two  possible  cases,  y  <  3  or  y  > 
3,  which  are  shown  in  Eigs.  4  and  5,  respeetively.  Because  y  =  1  +  2m/M,  note  that  y  <  3 
or  y  >  3  eorresponds  respeetively  to  simply  p  <  1  (m  <  M)  or  p  >  1  (m  >  M). 
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Figure  4.  Qualitative  sketch  (not  to  scale)  of  the  amplitude-frequency  plane  for  parametric  excitation  of  the 
lower-frequency  (^)  mode  by  the  upper-frequency  (p)  mode  with  amplitude  A.  This  case  corresponds  to 
m  <  M.  The  excitation,  and  thus  instability  of  the  p  mode,  occurs  for  points  on  the  heavy  line  that  he  in  the 

hatched  region. 


Figure  5.  Qualitative  sketch  (not  to  scale)  of  the  amplitude-frequency  plane  for  parametric  excitation  of  the 
lower-frequency  (^)  mode  by  the  upper-frequency  (p)  mode  with  amplitude  A.  This  case  corresponds  to 
m  >  M  (in  contrast  to  m  <  M  in  Fig.  3).  The  excitation,  and  thus  instability  of  the  p  mode,  is  not  predicted 
to  occur,  because  no  points  on  the  heavy  line  lie  in  the  hatched  region. 
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There  is  potentially  an  upper  threshold  in  both  oases.  However,  as  in  the  previous 
ease  (See.  II.B),  this  is  shown  below  to  oorrespond  to  O)  =  0,  and  so  is  predicted  to  have 
no  effect.  Due  to  the  relative  values  of  the  slopes  of  the  lines  of  the  V  and  the  amplitude- 
dependent  frequency,  a  hardening  (a  <  0)  upper-frequency  mode  with  amplitude  greater 
than  that  corresponding  to  the  dashed  line  in  Fig.  4  is  predicted  to  be  unstable,  which 
occurs  only  for  y  <  3.  For  y  >  3,  Fig.  5  shows  that  the  mode  is  not  subject  to  the 
instability. 

To  calculate  the  thresholds  of  the  instability,  we  substitute  the  amplitude- 
dependent  frequency,  Eq.  (12),  into  the  parametric  excitation  condition,  Eq.  (14),  and 
proceed  similarly  as  in  the  previous  case  (Sec.  II.B).  The  upper  threshold  value  of  the 
ordinate  in  Eigs.  4  and  5  is  then  found  to  be  aA^  =  dcOo^/Sy^,  which  corresponds  to  zero 
frequency  according  to  Eq.  (12).  The  theory  thus  predicts  that  this  threshold  does  not 
exist.  The  calculation  of  the  lower  threshold  in  Eig.  4  yields  the  condition 


-oA^ 


> 


4  Y-1  2 


3(1  +  2p)  (1-p) 


-(0„ 


(15) 


which  is  valid  for  y  <  3  or  p  <  1 .  The  theory  thus  predicts  that  a  hardening  (a  <  0)  upper- 
frequency  mode  is  unstable  for  amplitudes  that  satisfy  Eq.  (15),  and  that  a  softening  (a  > 
0)  upper-frequency  mode  is  not  subject  to  the  instability. 

As  in  the  previous  case,  the  results  are  based  only  on  a  leading-order  perturbation 
theory.  The  lines  in  Eigs.  4  and  5  are  actually  approximations  to  curves,  so  an  upper 
threshold  may  exist  for  O)  >  0.  We  resolve  this  with  numerical  simulations  in  Ch.  III. 
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III.  NUMERICAL  SIMULATIONS 


Numerical  simulations  of  the  exact  equations  of  motion,  Eqs.  (6)  and  (7),  were 
run  utilizing  Microsoft  Visual  Studio  1999  in  the  C++  programming  language.  To  view 
programs,  see  Appendix  A.  Both  the  upper-frequency  mode,  r),  and  lower-frequency 
mode,  were  examined  in  both  the  hardening  and  softening  nonlinear  cases. 
Dimensionless  units  were  used  where  (Oo  =  1  and  a  =  ±1.  All  simulations  were 
performed  using  the  Euler-Cromer  method.  (Garcia)  To  determine  the  stability  of  each 
mode  in  the  simulation  where  the  maximum  amplitude  was  unity,  the  mode  under 
investigation  was  given  an  initial  amplitude  and  the  other  mode  was  perturbed.  The 
perturbation  amplitude,  epsilon,  ranged  from  0.001  to  0.00001,  time  steps  ranged  from  5t 
=  0.01  to  0.0001  and  the  maximum  time  examined  ranged  from  tmax  =  100  to  5000.  Time 
series  were  imported  into  Microcal  Origin  6.0  and  plotted  for  observation.  Due  to  a 
decreasing  period  of  motion  for  larger  amplitudes  in  a  hardening  system,  a  sufficiently 
small  time  step  had  to  be  used  to  ensure  that  a  false  instability  was  not  documented.  An 
example  of  this  is  the  hardening  upper- frequency  mode  (Sec.  III.D)  for  p  =  m/M  >  1, 
which  we  originally  found  to  be  unstable  above  an  amplitude  threshold,  but  which  we 
later  realized  was  a  numerical  artifact  due  to  the  time  step  being  too  large.  The  criterion 
for  determining  stability  was  a  0.001  increase  in  the  amplitude  of  the  zero  amplitude 
mode.  In  agreement  with  the  theory,  only  the  hardening  upper- frequency  mode  and  the 
softening  lower-frequency  mode  displayed  the  parametric  instability. 

Additionally,  each  softening  mode  was  examined  for  the  point  at  which  the 
motion  would  diverge  if  the  system  were  provided  enough  energy  to  exhibit  “over-the- 
top”  motion.  This  can  be  viewed  similar  to  a  ball  in  a  valley  with  peaks  on  either  side.  If 
given  to  much  energy,  the  amplitude  of  oscillation  will  exceed  the  height  of  one  of  the 
peaks  and  the  ball  will  roll  over  the  peak  and  never  return.  We  call  this  the  “over-the- 
top”  amplitude,  Aott- 
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A,  SOFTENING  LOWER-FREQUENCY  MODE 

In  the  lower-frequency  mode,  the  case  where  the  coupling  softens  is  predicted  by 
the  theory  to  be  parametrically  unstable.  When  the  lower-frequency  mode  is  initiated 
with  great  enough  amplitude,  the  system  demonstrates  a  beating  motion  where  energy  is 
transferred  from  one  mode  to  the  other  (Fig.  6).  The  nonlinearity  causes  the  initially- 
excited  mode  to  parametrically  drive  the  other  mode  at  twice  the  frequency  of  the 
initially-excited  mode,  as  seen  in  Eq.  (9).  This  is  known  as  the  principal  parametric 
resonance,  where  the  frequency  of  the  response  locks  onto  one-half  the  frequency  of  the 
drive,  and  is  the  basis  of  our  theory.  In  the  general  theory  of  parametric  excitation,  see 
references  in  Ref.  (Denardo),  it  is  shown  that  there  occur  other  resonances  in  which  the 
ratio  of  the  response  to  drive  frequency  is  1  (isochronous  resonance),  3/2,  2,  5/2,  etc., 
although  these  are  successively  more  difficult  to  obtain. 
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Figure  6.  (a)  Unstable  oseillation  repeating  in  the  lower-frequency  mode  (gray  lines)  with  softening 

nonlinearity  for  p  =  0.1,  A  =  0.3267  and  criterion  =  0.001.  Space  between  oscillations  is  so  small  that 
oscillations  appear  as  gray  background.  Starting  amplitude  is  2%  greater  than  the  amplitude  threshold  for 
this  system,  (b)  Magnification  of  (a)  to  show  the  increasing  amplitude  in  the  upper-frequency  mode. 


When  the  initial  amplitude  of  the  lower-frequency  mode  is  decreased,  the 
resultant  amplitude  of  the  upper-frequency  mode  is  also  decreased.  Additionally,  the 
time  required  for  the  peak  to  occur  in  the  upper-frequency  mode  increases  (Fig.  7).  As 
the  amplitude  threshold  is  approached,  the  peak  amplitude  of  the  upper-frequency  mode 
approaches  zero  and  the  time  increases  to  infinity.  Thus,  it  may  appear  impossible  to 
accurately  determine  the  amplitude  threshold.  However,  sufficient  accuracy  on  the  scale 
of  Fig.  8  (see  below)  is  not  difficult  to  obtain. 
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Figure  7.  Unstable  oscillation  in  the  lower-frequency  mode  (gray  lines)  with  softening  nonlinearity  for  p  = 
0.1,  A  =  0.3204  and  criterion  =  0.001.  Space  between  oscillations  is  so  small  that  oscillations  appear  as 
gray  background.  The  threshold  amplitude  for  this  system  is  Ath  =  0.3203. 

The  dependence  between  the  threshold  amplitude  for  the  parametric  instability 
and  the  mass  ratio  of  the  system  can  be  seen  in  Fig.  8.  As  the  mass  ratio  increases,  so 
does  the  amplitude  threshold.  The  theory  is  only  expected  to  be  valid  when  the 
frequencies  of  the  normal  modes  are  near  each  other,  which  corresponds  to  p  «  1 .  The 
circles  correspond  to  the  amplitude  threshold  found  through  numerical  simulation  for 
values  of  p  ranging  from  0.01  to  2.  These  values  are  plotted  against  equation  (11)  and 
are  in  good  agreement  for  p  <  0.98.  However,  at  p  □  0.98  the  system  no  longer  appears 
to  exhibit  the  instability  but  goes  straight  into  “over-the-top”  motion  given  great  enough 
amplitude.  This  is  in  sharp  contrast  to  the  theoretical  over-the-top  amplitude  of  1 . 

Figure  8  also  shows  the  amplitude  threshold  for  divergent  (unbounded)  motion. 
For  amplitudes  between  the  circles  and  the  squares  the  parametric  instability  leads  to  the 
increasing  beating  amplitudes  and  increasing  beating  frequency  described  above.  The 
divergent  motion  can  occur  because  the  coupling  force,  Eq.  (1),  for  a  >  0  corresponds  to 
a  potential  energy  that  is  a  well  only  for  |x|  <  (k/am)  =  (OoOC  ,  and  negatively  diverges 
as  |x|  ^  00.  The  central  mass  M  is  at  rest  for  the  lower-frequency  normal  mode,  so  the 
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1/2 

expected  amplitude  threshold  for  this  “over-the-top”  instability  is  simply  cooOC  =  1  (in 
dimensionless  units).  The  data  in  Fig.  8  show  that  the  parametric  instability  causes  over- 
the-top  motion  to  occur  for  amplitudes  that  are  substantially  less  than  the  expected  value, 
unless  p  «  1.  This  behavior  is  not  yet  understood,  and  is  in  contrast  to  the  case  of 
linearly  coupled  nonlinear  oscillators.  (Denardo)  For  p  ^  1,  the  parametric  instability  at 
the  amplitude  threshold  initiates  the  over-the-top  motion,  to  at  least  three  significant 
figures  in  the  value  of  the  amplitude. 


Ratio  of  Masses,  p  =  m/M 


Figure  8.  Deviation  between  the  theoretical  calculation  and  the  numerical  simulation  results  for  the 
threshold  amplitude  of  the  lower-frequency  mode  with  softening  nonlinearity.  Additionally,  the  “over-the- 
top”  amplitude  approaches  the  threshold  amplitude  and  dominates  at  p  >  1.0.  Ajj,  =  Athreshoid 


Data  was  gathered  for  values  of  p  beyond  those  in  Fig.  8,  and  discontinuities  in 
the  threshold  amplitude  as  a  function  of  p  were  found  (Fig.  9).  The  first  discontinuity  is 
at  p  =  3.0,  where  the  threshold  amplitude  jumps  from  0.71  to  0.53  as  p  is  increased.  The 
second  is  at  p  =  7.1,  where  the  jump  is  from  0.64  to  0.54.  Examination  of  time  series 
reveals  that  the  first  discontinuity  corresponds  to  the  isochronous  parametric  resonance. 
This  continues  to  occur  until  the  next  discontinuity,  which  corresponds  to  the  3/2 
parametric  resonance.  For  this  range  of  p  values,  the  principal  parametric  resonance 
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presumably  continues  to  exist,  but  the  threshold  for  the  isochronous  resonance  is  smaller. 
There  is  apparently  an  infinite  sequence  of  discontinuities  that  corresponds  to  the 
sequence  of  parametric  instabilities  described  above.  For  examples  of  the  sequence  of 
parametric  instabilities  see  Appendix  B. 


Figure  9.  Instabilities  due  to  higher  order  parametric  resonances  in  the  lower-frequency  softening  mode  at 

high  values  of  p.  A,h  =  Athreshoid 

In  Sec.  II. B,  the  possibility  was  mentioned  of  an  upper  threshold  existing  for  O)  > 
0  due  to  the  lines  in  Fig.  3  being  approximations  to  curves.  However,  numerical 
simulation  showed  no  such  threshold  for  0  <p<  1.  The  instability  may  be  prevented  by 
the  presence  of  the  over-the-top  instability  of  our  model  force,  Eq.  (1).  Other  softening 
forces  that  do  not  possess  this  property  may  yield  an  upper  threshold.  An  example  is  the 
force  whose  potential  energy  in  dimensionless  units  is  ln(l+x  ). 

B,  HARDENING  LOWER-FREQUENCY  MODE 

The  theory  also  predicts  that  the  hardening  lower-frequency  mode  should  be 
stable  for  the  parametric  instability.  Numerical  simulations  show  that  this  mode  is  stable 
for  p  <  0.5,  except  at  p  ~  0  due  to  the  curve  approximation  of  Fig.  3.  For  0.5<p<  1.2, 
the  system  enters  a  region  of  instability  which  examination  of  the  time  series  shows  is 
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caused  by  the  isoehronous  parametrie  resonanee,  App.  B.  The  instability  repeats  for  2.0 
<  p  <  3.1  and  4.1  <  p  <  5.0.  Examination  of  the  time  series  showed  that  these  regions 
were  the  and  4*'’  order  parametrie  resonanees,  respeetively.  Thus,  we  were  able  to 
isolate  the  regions  of  E*  through  4*’’  order  parametrie  resonance  between  0  <  p  <  5,  see 
Fig.  10.  The  figure  elearly  shows  regions  of  stability  between  regions  of  instability.  It 
appears  that  as  p  is  inereased,  sueeessive  higher-order  parametrie  resonanees  oeeur  over 
intervals  that  are  separated  by  intervals  of  stability.  Care  had  to  be  taken  in  ensuring  the 
time  steps  were  low  enough  that  no  signifieant  ehange  oeeurred  in  the  motion  between 
time  steps.  It  should  be  noted  that  “over-the-top”  motion  is  not  possible  for  a  hardening 
nonlinearity. 


Figure  10.  Mass  dependence  on  the  amplitude  threshold  of  the  lower-  frequency  mode  with  hardening 

nonlinearity. 

As  mentioned  in  See.  (II. B),  the  possibility  exists  for  a  region  of  instability  near  p 
~  0,  beeause  the  lines  in  Fig.  3  are  actually  approximations  to  curves.  Numerical 
simulations  show  that  for  0.001  <p<0  .009  the  instability  does  exist,  see  Figs.  10  and  11. 
The  amplitude  threshold  increases  monotonically  from  zero  as  p  inereases,  and  has  the 
value  1.1591  at  p  =  0.005. 
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Figure  11.  Threshold  amplitude  for  the  lower-  frequency  mode  with  hardening  nonlinearity  at  low  p  values. 


C.  SOFTENING  UPPER-FREQUENCY  MODE 

According  to  the  theory  in  See.  II. C,  the  upper-frequeney  mode  with  softening 
nonlinearity  was  predieted  by  theory  not  to  be  subjeet  to  the  parametrie  instability. 
Additionally,  over-the-top  motion  should  oeeur  when  the  ehange  in  length  of  a  spring 
exeeeds  the  value  eorresponding  to  zero  foree  in  Eq.  (1),  whieh  is  (Oo/oc  .  For  the  upper- 
frequeney  mode  with  amplitude  A,  the  eentral  mass  M  moves  in  antiphase  with  amplitude 
2pA,  so  the  ehange  in  length  of  a  spring  is  (1  +  2p)A  =  yA.  The  theoretieal  maximum 
stable  amplitude  of  the  upper  mode  is  thus  rimax  =  coo/yot  =  l/y  (in  dimensionless  units). 

However,  numerieal  simulations  revealed  growth  in  the  lower  mode  for  0  <  p  < 
0.017,  and  over-the-top  motion  for  all  values  of  p  (Fig.  12).  However  the  “over-the-top” 
amplitude  was  10-20%  less  than  the  theoretieal  value,  depending  upon  the  value  of  y 
unless  y  ~  1  (p  ~  0).  The  diserepaney  persists  even  when  we  make  the  amplitude  of  the 
perturbation  of  the  lower  mode  extremely  small. 
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Figure  12.  Amplitude  threshold  for  the  upper-frequency  mode  with  softening  nonlinearity  Ath  =  Athreshoid- 

Zooming  into  the  point  of  initial  perturbation  shows  isoehronous  resonanee  and 
eontinuous  inerease  in  the  amplitude  of  the  lower-frequeney  mode  until  the  system 
exhibits  over-the-top  motion  (Fig.  13).  When  this  happens,  the  upper- frequeney  mode 
exhibits  over-the-top  motion  as  well.  The  direetion  of  inereasing  amplitude,  positive  or 
negative,  depends  on  the  direetion  of  perturbation.  The  interesting  point  to  note  is  that 
the  amplitude  of  the  lower  mode  never  goes  below  the  value  of  the  perturbation  it 
reeeived,  it  either  stays  at  values  greater  than  0.0001  or  less  than  -0.0001  (Fig.  14). 
Thus,  it  appears  that  the  isoehronous  parametrie  resonanee  is  responsible  for  initiating  the 
over-the-top  motion.  This  instability  of  the  softening  upper  mode  is  in  eontrast  to  the 
previously  investigated  ease  of  two  linearly  eoupled  nonlinear  oseillators,  (Denardo)  for 
whieh  it  was  found  that  the  mode  was  stable  for  all  amplitudes  less  than  the  theoretieal 
over-the-top  amplitude. 
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Figure  13.  “Over-the  top”  motion  in  the  upper- frequency  mode  with  softening  nonlinearity. 
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Figure  14.  Zoom  of  “over-the-top”  motion  in  the  upper-frequency  mode  with  softening  nonlinearity 


As  explained  in  See.-  II. C,  when  p  ~  0  there  is  a  possibility  of  a  threshold 
amplitude  above  whieh  the  softening  upper- frequency  mode  may  be  subject  to  the 
principal  parametric  instability.  Numerical  simulations  show  that  this  does  occur  for 


23 


0<p<0  .016  (Fig.  15).  The  amplitude  threshold  inereases  monotonically  from  zero  as  p 
increases,  and  has  the  value  0.662  at  p  =  0.005.  The  threshold  smoothly  connects  with 
the  isochronous  resonance  threshold  at  p  =  0.017.  For  p  not  near  zero,  the  existence  of 
principal  parametric  resonance  cannot  be  resolved  due  to  the  presence  of  the  over-the-top 
instability  of  our  force  model,  Eq.  (1),  as  in  the  case  of  the  softening  lower- frequency 
mode  (Sec-  UFA) 
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Figure  15.  Amplitude  threshold  for  the  upper- frequency  mode  with  softening  nonlinearity  in  the  region  of 

p  ~  0.  Atji  —  Athreshold- 

D.  HARDENING  UPPER-FREQUENCY  MODE 

In  the  case  where  the  coupling  hardens,  the  theory  in  Sec.  11. C  predicts  the  upper- 
frequency  mode  to  be  unstable  for  the  parametric  instability  for  amplitudes  above  the 
threshold  value  given  by  Eq.  (15)  for  0  <  p  <  1,  and  to  not  exhibit  the  instability  for  p  >  1. 
As  with  the  lower-frequency  mode  with  softening  nonlinearity,  when  the  upper- 
frequency  mode  is  initiated  with  great  enough  amplitude,  the  system  demonstrates  the 
same  beating  motion  where  the  frequency  of  the  two  modes  is  the  same  and  is  consistent 
with  the  principal  parametric  resonance  (Eig.  16).  The  greater  the  amplitude  above 
threshold,  the  higher  the  repetition  frequency. 
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Figure  16.  Unstable  oscillation  repeating  in  the  upper- frequency  mode  (gray  lines)  with  hardening 
nonlinearity  for  p  =  0.1,  A  =  0.3241  and  criterion  =  0.001.  Space  between  oscillations  is  so  small  that 
oscillations  appear  as  gray  background.  Starting  amplitude  is  1%  greater  than  the  amplitude  threshold  for 

this  system. 

Numerical  simulations  of  the  exact  equations  of  motion,  Eqs.  (6)  and  (7), 
identified  the  instability  with  increasing  amplitude  of  the  lower-frequency  mode  from 
zero,  as  seen  in  Fig.  17.  When  the  initial  amplitude  of  the  upper- frequency  mode  is 
decreased,  the  resultant  amplitude  of  the  lower-frequency  mode  is  also  decreased. 
Additionally,  the  time  required  for  the  peak  to  occur  in  the  lower-frequency  mode  also 
increases.  As  the  amplitude  threshold  is  approached,  the  peak  amplitude  of  the  upper- 
frequency  mode  approaches  zero  and  the  time  increases  to  infinity. 
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Figure  17.  Unstable  oscillation  in  the  upper- frequency  mode  (gray  lines)  with  hardening  nonlinearity  for  p  = 
0.1,  A  =  0.3209  and  criterion  =  0.01.  Space  between  oscillations  is  so  small  that  oscillations  appear  as  gray 
background.  The  threshold  amplitude  for  this  system  is  0.3208. 


The  dependence  between  the  threshold  amplitude  and  the  mass  ratio  can  be  seen 
in  Fig.  18.  As  the  mass  ratio  increases,  so  does  the  amplitude  threshold.  The  amplitude 
threshold  was  found  through  numerical  simulations  for  0  <  p  <  1.  The  system  did  not 
exhibit  the  instability  for  1  <  p  <  2.  The  experimental  values  were  plotted  against  Eq. 
(15)  and  are  in  good  agreement  for  p  <  1,  even  at  p  =  0.999  where  the  approximate 
analytical  amplitude  threshold  is  only  1.3%  less  than  the  numerical  simulation  value. 


26 


Figure  18.  Theoretical  calculation  and  the  numerical  simulation  results  for  the  threshold  amplitude  of  the 
upper-frequency  mode  with  hardening  nonlinearity.  Ath  =  Athreshoid- 

Based  on  our  experience  with  the  previous  mode  (Sec.-  III.A),  we  expected  the 
parametric  resonances  other  than  the  principal  resonance  would  occur  for  p  >  1 .  This  is 
apparently  not  the  case,  however,  because  numerical  simulations  for  a  wide  variety  of 
mass  ratio  and  amplitude  values  yield  a  stable  upper-frequency  mode,  as  long  as  the  time 
step  is  chosen  to  be  sufficiently  small. 

As  explained  in  Sec.-  II. C,  the  hardening  upper-frequency  mode  for  p  ~  0  may 
become  stable  above  an  upper  amplitude  threshold.  However,  numerical  simulations 
show  no  evidence  of  this. 
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IV.  DEMONSTRATIONS 


A,  THREE  COUPLED  PENDULUM  SYSTEM 

For  a  visual  example  of  the  two  linear  normal  oseillatory  modes  and  the 
instability  of  the  triatomie  moleeule,  the  triatom  system  was  translated  into  a  three- 
eoupled  pendulum  system,  see  Fig.  19.  This  system  is  similar  to  the  two-eoupled- 
pendulum  system  apparatus  diseussed  in  Ref.  (Denardo).  This  system  has  three  linear 
normal  modes.  In  the  fundamental  normal  mode,  all  three  pendulums  oseillate  in  phase 
with  the  same  amplitude.  The  two  other  linear  normal  modes  are  similar  to  the  two 
oseillatory  normal  modes  of  the  triatomie  moleeule,  the  lower  and  upper  modes.  There  is 
a  distinetion  between  the  fundamental  mode  of  the  pendulum  system  and  the  zero- 
frequeney  mode  of  the  triatomie  moleeule,  but  this  is  irrelevant  in  regard  to  the 
oseillatory  modes  of  the  triatomie  moleeule.  There  is  also  a  distinetion  in  the 
nonlinearity  of  the  two  systems.  The  pendulums  are  eoupled  nonlinear  oseillators, 
whereas  the  triatomie  moleeule  is  a  system  of  nonlinearly  eoupled  masses. 


Figure  19.  Three-coupled  pendulum  system  where  the  center  mass  is  heavier  than  the  outer  masses. 
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Pendulums  #1  &  #3  were  of  equal  size  and  weight,  m  =  39.5  g.  For  pendulum  #2, 
aluminum,  brass  and  lead  weights  were  used.  Pendulum  #2  was  mAi  =  79.3  g,  Pai  = 
0.498,  mb  =  246.2  g,  Pb  =  0.161,  and  mpb  =  302.1  g,  ppb  =  0.131.  The  distanee  from  the 
eenter  of  the  large  V  is  di  =  12  em.  The  distanee  from  the  eenter  of  the  small  V  is  d2  =  22 
em.  The  length  of  the  large  V  is  h  =  39.5  em.  The  length  of  the  small  V  is  b  =  10  em. 
The  aluminum  weight  was  employed  in  the  lower  and  upper  frequeney  modes  at  low  and 
high  amplitudes  and  all  four  variations  of  the  system  were  stable.  The  same  was  true  for 
the  brass  and  lead  weight  system.  This  is  perhaps  not  surprising  beeause  the  nonlinearity 
of  the  two  systems  is  different  (see  above).  However,  we  did  find  that  the  fundamental 
mode  of  the  three-pendulum  system  is  unstable  above  an  amplitude  threshold.  This  is  the 
few-degree-of-freedom  limit  of  the  Benjamin-Feir  instability  in  eontinuous  systems. 
(Denardo) 

B,  MAGNASWING 

Next,  due  to  its  similarity  to  the  triatomie  moleeule,  the  eommereial 
MagnaSwing^"^  toy  (Denardo)  was  investigated  using  three  magnetie  pendulums  set  to 
repel  one  another  (Fig.  20).  As  sold,  this  is  a  system  of  as  many  as  five  pendulums  that 
oseillate  longitudinally,  where  magnets  in  the  bobs  supply  the  eoupling.  As  in  the  above 
system  of  eoupled  pendulums,  the  fundamental  mode  is  distinet  from  the  zero-frequeney 
mode  of  the  triatomie  moleeule,  but  this  is  irrelevant.  This  system  does  display  two  other 
modes  whieh  are  similar  to  the  upper  and  lower  frequeney  modes  of  the  triatomie 
moleeule  system. 

Additionally,  the  nonlinearity  is  distinet  from  the  triatomie  moleeule  not  only  as  a 
result  of  the  pendulum  nonlinearity,  but  also  beeause  the  magnetie  eoupling  is  modeled 
by  both  eubie  and  quadratie  nonlinearities  in  the  foree.  Our  theory  (Sees.-  II. B  and  II. C) 
assumes  a  pure  eubie  nonlinearity.  It  is  thus  perhaps  not  surprising  that  the  two  modes 
that  are  similar  to  the  triatomie  oseillatory  modes  did  not  exhibit  the  instability.  The 
pendulums  were  all  of  equal  size  (d=  32  mm)  and  weight  (m  =  41.9  g).  The  eenter 
pendulum  was  weighted  with  a  solid  brass  eylinder  (d  =  30.75  em,  m  =  106.1  g)  to 
aehieve  a  mass  ratio  of  p  =  0.284.  As  in  the  eoupled  pendulum  system,  even  with  a 
substantially  inereased  eentral  mass,  the  instability  still  did  not  oeeur.  The  lower  mode  of 
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vibration  was  stable  at  all  amplitudes.  The  upper  mode  of  vibration,  though  diffieult  to 
find,  was  also  stable  at  all  amplitudes. 


Figure  20.  A  three-pendulum  configuration  of  the  MagnaSwing  toy.  The  pendulums  were  all  the  same  size 
and  weight  and  additional  solid  brass  weight  was  added  to  the  center  pendulum  to  achieve  a  mass  ratio  of 

0.284. 

C.  HARDENING  AIR  TRACK  GLIDER  SYSTEM  (MODEL  1) 

The  third  system  I  investigated  was  three  coupled  gliders  on  an  air  track  (Fig.  21). 
The  gliders  were  connected  by  two  springs  (one  on  each  side  to  maintain  symmetry)  and 
mounted  with  repulsive  ceramic  magnets.  The  purpose  of  the  magnets  was  to  introduce  a 
hardening  nonlinearity.  This  system  is  identical  to  the  triatomic  model  with  the  exception 
that  the  coupling  can  be  approximated  by  both  quadratic  and  cubic  nonlinearities, 
whereas  the  model  assumes  a  pure  cubic  nonlinearity.  The  mass  of  the  magnets  was 
between  47.7  g  and  49.0  g.  Each  glider  was  weighed  and  multiple  combinations  of 
magnets  were  used  to  ensure  gliders  #1  &  #3  were  always  of  the  same  weight.  Multiple 
variations  of  this  system  were  examined  using  gliders  with  mass  150  g,  300  g  and  400  g. 
Mass  ratios  between  0.269  and  0.8  were  achieved.  Neither  of  the  oscillatory  normal 
modes  were  found  to  exhibit  the  instability,  even  at  large  amplitudes  and  small  mass 
ratios. 
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Figure  21.  Air  track  glider  setup  to  mimic  the  longitudinal  vibration  of  a  triatomic  molecule  using  three 

spring  coupled  gliders  and  ceramic  magnets. 
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In  the  demonstrations  deseribed  so  far,  it  should  be  noted  that  the  lower-frequency 
mode  can  be  readily  excited  by  releasing  the  outer  masses  simultaneously  from  rest  at  the 
same  amplitude  in  opposite  directions,  with  the  central  mass  at  rest.  In  the  first  two 
demonstrations,  pure  upper-frequency  mode  motion  can  be  difficult  to  excite  by 
attempting  to  release  the  masses  with  the  proper  relative  amplitudes.  An  alternative  is  to 
use  a  finger  to  drive  a  mass  or,  in  the  case  of  a  pendulum,  a  point  near  the  pivot.  By 
resonantly  driving  one  of  the  masses  at  the  frequency  of  the  upper  mode,  patiently 
waiting  for  transient  motion  to  sufficiently  decay,  and  carefully  removing  the  drive,  an 
approximately  pure  upper  mode  can  be  obtained.  In  the  third  demonstration,  as  well  as 
the  demonstrations  below,  the  upper-frequency  mode  is  easily  excited  by  simply  holding 
M  at  rest  and  displacing  the  two  m’s  the  same  distance  in  the  same  direction,  and 
simultaneously  releasing  all  the  masses  from  rest.  As  long  as  the  system  is  symmetric, 
this  initial  condition  excites  only  the  upper-frequency  mode. 

Although  none  of  the  above  four  systems  exhibit  the  desired  instability,  it  should 
be  noted  that  each  of  these  apparatus  offers  an  easy  means  of  demonstrating  the  two 
linear  normal  oscillatory  modes  of  the  triatomic  molecule.  This  is  useful  in  the  classroom 
because  it  offers  a  concrete  visualization  of  the  motion. 

D,  HARDENING  AIR  TRACK  GLIDER  SYSTEM  (MODEL  2) 

Our  failure  with  the  quick  demonstrations  described  above  prompted  us  to  build 
apparatus  with  symmetric  potential  wells,  for  which  the  theory  in  Secs.-  IIB  and  II. C  is  an 
approximation  due  to  the  cubic  correction  to  a  linear  force.  In  contrast  to  the  simple 
apparatus  for  the  previously-investigated  instability  of  two  linearly  coupled  nonlinear 
oscillators,  (Denardo)  the  apparatus  here  require  much  more  effort  to  construct.  Figure 
22  shows  a  system  with  an  air  track,  gliders,  springs  and  magnets  that  are  realizations  of  a 
classical  triatomic  molecule  with  symmetric  nonlinear  potentials.  The  nonlinearity  of  the 
apparatus  is  due  to  the  repulsive  forces  of  the  magnets,  and  the  outer  magnets  on  the 
outer  gliders  are  for  balance  only,  in  order  to  minimize  friction.  We  use  strong  ceramic 
magnets  of  diameter  2.2-cm  and  length  2.5-cm.  The  magnets  are  taped  to  7.5-cm  long 
channels  made  by  cutting  PVC  pipe  lengthwise  in  half.  The  channels  are  attached  to  the 

gliders  with  screws.  In  this  system  the  upper-frequency  mode  should  be  unstable  as  long 
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as  p  =  m/M  <  1.  The  outer-most  gliders  were  added  to  the  ends  of  model  1  and 
conneeted  by  a  1.9  em  diameter  and  152.4  em  long  thin  aluminum  tube  to  ereate  a  single 
mass  M,  formed  by  the  middle  and  outer  gliders.  This  pulled  gliders  #2  and  #4  further 
away  from  the  eenter  glider.  Additionally,  beeause  gliders  #1,  #3  and  #5  were  rigidly 
attached,  the  mass  of  the  three  gliders  combined  to  give  the  larger  mass,  M,  and  enabled 
the  system  to  achieve  lower  mass  ratios. 
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Figure  22.  Double-spring  air  track  glider  system  setup  to  mimic  a  triatomic  molecule  using  five  spring 
coupled  gliders.  Gliders  #1,  #3,  #5  are  rigidly  coupled  to  form  mass  M. 


The  tube  is  attached  to  the  gliders  by  posts  with  washers,  where  it  is  important  to 
note  that  the  washers  must  be  carefully  inserted  or  removed  such  that  the  entire  M  mass 
moves  with  little  friction.  Without  such  an  adjustment,  an  excessive  downward  force  due 
to  the  aluminum  tube  can  cause  a  glider  to  stick  on  the  track.  In  addition,  the  attached 
gliders  must  be  free  to  swivel,  so  that  friction  is  minimized.  We  use  a  Daedalon  2.5  m  air 
track  with  a  gold-colored  (0.15  kg)  for  each  m  mass  and  three  blue  (0.45  kg)  gliders 
connected  together  for  the  M  mass.  We  also  use  Daedalon  2.5  cm  springs.  (Daedalon) 
The  air  track  is  driven  by  two  air  pumps,  one  attached  to  each  end. 

Another  particular  problem  involves  symmetry.  If  a  system  is  not  symmetric, 
then  releasing  the  gliders  in  the  simple  manner  described  above  will  not  lead  to  pure 
excitation  of  a  normal  mode,  so  linear  beating  will  occur  between  the  normal  modes. 
This  is  unacceptable  in  our  case,  because  the  signature  of  the  parametric  instability  is  a 
nonlinear  beating  between  the  modes.  Due  primarily  to  variations  in  the  spring  constants 
of  the  springs,  we  expended  a  substantial  amount  of  effort  in  arranging  the  system  in  Fig. 
22  to  be  sufficiently  symmetric.  Our  test  for  the  symmetry  was  the  absence  of  linear 
beating  for  our  simple  initial  conditions. 
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In  the  initial  investigation  the  gliders  remained  eonneeted  by  two  springs,  one 
mounted  on  each  side  of  all  five  gliders.  Mass  ratios  between  0.269  and  0.8  were 
achieved.  It  was  found  that  neither  of  the  oscillatory  normal  modes  exhibited  the 
instability,  even  at  large  amplitudes  and  small  mass  ratios. 

In  the  second  system  each  of  the  gliders  were  coupled  by  a  single  2.5  cm  spring. 
Figure  23  is  a  picture  of  the  actual  apparatus  developed.  This  was  a  revised  version  of 
the  double-spring  air  track  glider  system  described  previously.  Holes  were  drilled  in  the 
gliders  to  attach  a  single  spring  2-mm  above  the  air  track  and  2-cm  below  the  magnets. 
The  separation  distance  between  gliders  was  21.25-cm.  At  large  amplitudes  of  the 
motion,  the  m  gliders  can  depart  from  the  air  track.  Because  the  threshold  amplitude  for 
the  instability  decreases  as  p  decreases,  the  value  of  p  must  be  sufficiently  small  such  that 
the  instability  occurs  before  the  gliders  depart  from  the  air  track,  but  not  so  small  that  the 
linear  motion  below  the  amplitude  threshold  cannot  be  clearly  demonstrated.  To  achieve 
this,  as  stated  earlier,  we  use  blue  gliders  for  the  M  mass  and  gold  gliders  for  the  m 
masses.  In  addition,  we  add  a  1 00-gram  mass  to  each  of  the  two  sides  of  each  outer 
glider.  The  resultant  mass  ratio  (which  includes  the  aluminum  tube,  posts,  and  magnets) 
is  p  =  0.1 1.  The  lower- frequency  mode  is  observed  to  be  stable  at  all  amplitudes,  where 
as  the  upper-frequency  mode  is  observed  to  be  stable  at  smaller  amplitudes  and  unstable 
at  larger  amplitudes. 


Figure  23.  Picture  of  actual  air  track  glider  system  utilizing  repulsive  magnets  and  coupled  with  single 
springs  to  demonstrate  the  instability  in  the  upper- frequency  hardening  mode,  p  =  0.1 1. 
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First,  gliders  #2  &  #4  were  displaced  2.0-cm  in  the  same  direction  followed  by  2- 
cm  in  the  opposite  direction  to  ensure  the  upper  and  lower  frequency  modes  were  stable 
at  low  amplitudes  and  the  system  is  symmetric.  Then,  gliders  #2  and  #4  were  displaced 
7.0  cm  in  opposite  directions  to  see  if  there  was  any  instability  in  the  lower  frequency 
mode,  there  was  none.  This  displacement  distance  was  chosen  because  it  was  the 
maximum  displacement  which  did  not  cause  chaotic  behavior  in  the  gliders,  resulting  in 
rubbing  between  the  glider  and  the  air  track  or  resulting  in  gliders  #2  &  #4  jumping  off 
the  air  track.  Thorough  examination  demonstrated  that  this  mode  of  vibration  was  stable 
even  when  the  gliders  were  not  displaced  exactly  equal.  Instead  of  this  mode  becoming 
unstable,  the  entire  system  would  also  exhibit  the  zero-frequency  mode  and  translate 
down  the  track.  Finally,  gliders  #2  and  #4  were  displaced  in  the  same  direction  7.0  cm  to 
look  for  the  instability.  After  thorough  examination,  p  =  0.127  was  discovered  to  be  the 
minimum  mass  ratio  to  exhibit  the  principal  parametric  instability  under  the  above 
conditions. 

E.  SOFTENING  (SLANTED-SPRING)  AIR  TRACK  GLIDER  SYSTEM 

Figure  24  shows  the  top  view  of  a  softening  system,  for  which  the  lower- 
frequency  mode  should  be  unstable  above  an  amplitude  threshold.  To  leading  nonlinear 
order,  a  calculation  shows  that  the  slanted  springs  give  rise  to  a  softening  cubic 
nonlinearity  in  the  force,  as  long  as  the  springs  are  sufficiently  nonperpendicular  to  the 
track  (App.  C).  Specifically,  the  distance  from  the  center  of  an  m  glider  to  the  attachment 
point  of  the  spring  on  the  transverse  rod  must  be  less  than  half  the  distance  between  the 
centers  of  an  m  glider  and  the  central  M  glider  in  equilibrium.  In  our  apparatus,  the 
springs  were  initially  at  roughly  45°.  We  did  not  observe  the  predicted  instability,  even 
for  the  maximum  amplitude  in  which  the  gliders  touch. 


35 


Calculation  showed  that  in  order  to  maximize  the  softening  force  in  the  system 
the  length,  b,  from  the  center  of  glider  m  to  the  attachment  point  of  the  spring  on  the 
transverse  rail  had  to  be  approximately  0.7a,  where  a  here  is  the  distance  from  the  end  of 
glider  M  to  the  center  of  glider  m  and  the  length  of  the  unstreached  spring,  c,  had  to  be 
approximately  1.2a,  see  Fig.  25.  In  this  figure,  the  ordinate  is  the  ratio  of  the  nonlinear 
coefficient  to  the  linear  coefficient  in  the  force.  (Refer  to  App.  C.)  Unfortunately,  the 
optimized  system  was  not  achievable  due  to  a  large  variability  in  the  Daedelon  springs, 
which  prevented  symmetrization.  This  was  true  even  though  the  springs  had  the  same 
length  and  radius,  were  made  with  the  same  diameter  of  wire,  and  came  from  the  same  lot 
number.  Pairs  of  2.5  and  7.5  cm  springs  were  used  for  the  system  in  Fig.  24. 
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Figure  25.  Calculation  of  maximum  softening  force  for  varying  lengths  in  the  slanted-spring  air  glider 

demonstration. 

The  lowest  value  of  b  obtained  for  this  system  while  maintaining  symmetry  in  the 
springs  was  b  =  16  cm.  The  threshold  amplitude  required  to  initiate  the  instability  in  this 
system  was  calculated  from  Eq.  (11),  where  (Oo  =  27tfo,  p  =  m/M,  and  a  =  (2k/m)(4a  - 

bV(a^+b^)^%^c,  a  =  23  cm,  b  =  16  cm,  c  =  16.8  cm,  m  =  178.5  g,  M  =  2624  g.  The 
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spring  constants  were  eomputed  to  be  k  =  951  dynes/em.  From  this  the  threshold 
amplitude  was  determined  to  be  14.5  cm  and  the  maximum  amplitude  aehievable  in  the 
system  was  18.5  em.  It  was  apparent  that  the  system  was  above  threshold,  however,  after 
two  oscillations  the  amplitude  deereased  below  the  threshold  amplitude.  This  may  not  be 
enough  time  for  the  instability  to  develop. 

To  increase  the  softening  nonlinearity  and  assist  the  system  in  remaining  above 
the  threshold  amplitude  for  indueing  the  instability,  attractive  magnets  were  added  to  the 
system  (Fig.  26).  The  same  magnet  assembly  used  in  See.  IV.E,  was  used  in  the  slanted- 
spring  assembly.  However,  to  keep  the  maximum  amplitude  of  the  system  at  the  same 
value  as  in  See.  IV.F  the  magnets  on  gliders  #1,  #3,  and  #5  were  recessed  toward  the 
center  of  each  glider  8.5  em.  Beeause  the  threshold  amplitude  for  the  instability 
decreases  as  p  deereases,  the  value  of  p  must  be  suffieiently  small  sueh  that  the  instability 
oecurs  before  the  gliders  hit,  but  not  so  small  that  the  linear  motion  below  the  amplitude 
threshold  eannot  be  elearly  demonstrated.  To  aehieve  this,  we  use  blue  gliders  for  the  M 
masses  and  gold  gliders  for  the  m  masses.  In  addition,  I  added  150  grams  of  mass  to  eaeh 
of  the  two  sides  of  each  outer  glider  and  200  grams  of  mass  to  each  of  the  two  sides  of 
the  middle  glider.  The  resultant  mass  ratio  (whieh  includes  the  aluminum  tube,  posts, 
and  magnets)  is  p  =  0.095.  The  lower- frequency  mode  is  observed  to  be  stable  at  all 
amplitudes,  wheras  the  upper-frequency  mode  is  observed  to  be  stable  at  smaller 
amplitudes  and  unstable  at  larger  amplitudes. 


Figure  26.  Picture  of  actual  air  track  glider  system  utilizing  attractive  magnets  and  coupled  with  two  slanted 
springs  to  demonstrate  the  instability  in  the  lower-frequency  hardening  mode,  p  =  0.095. 
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First,  gliders  #2  &  #4  were  displaeed  2.0  em  in  the  same  direetion  followed  by  2 
em  in  the  opposite  direetion  to  ensure  the  upper  and  lower  frequeney  modes  were  stable 
at  low  amplitudes  and  the  system  is  symmetrie.  Then,  gliders  #2  and  #4  were  displaeed 
18.0  em  in  the  same  direetion  to  eonfirm  that  there  was  no  instability  in  the  upper- 
frequeney  mode.  This  displaeement  distanee  was  ehosen  beeause  it  was  the  maximum 
displaeement  that  did  not  eause  eollisions  between  gliders,  resulting  in  deereased 
amplitudes  between  gliders  and  indueed  beating.  Thorough  examination  demonstrated 
that  this  mode  of  vibration  was  stable  even  when  the  gliders  were  not  displaeed  exaetly 
equal.  Instead  of  this  mode  beeoming  unstable,  the  entire  system  would  also  exhibit  the 
zero  frequeney  mode  and  translate  down  the  traek.  Finally,  gliders  #2  and  #4  were 
displaeed  in  opposite  direetions  18.0  em  to  look  for  the  instability.  The  prineipal 
parametrie  instability  appeared  evident  but  was  very  weak.  To  inerease  the  instability  on 
the  system,  the  mass  ratio  was  redueed  to  p  =  0.089  by  adding  an  additional  200  grams  of 
mass  to  eaeh  of  the  two  sides  of  the  end  and  middle  gliders.  Contrary  to  the  initial 
assumption,  the  instability  deereased.  The  reason  for  this  behavior  is  that  deereasing  p 
deereases  the  threshold  amplitude,  but  longer  time  is  required  for  the  instability  to 
develop. 
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V.  CONCLUSIONS 


Due  to  nonlinearity,  a  normal  oscillatory  mode  of  a  system  can  parametrically 
drive  another  normal  oscillatory  mode.  Through  conservation  of  energy,  the  initially 
excited  mode  of  a  nonlinearly  coupled  system  must  be  unstable.  We  have  investigated  a 
parametric  instability  of  a  simple  model  of  a  classical  triatomic  molecule.  The  system  is 
rectilinear  and  symmetric,  and  the  coupling  includes  a  cubic  nonlinearity  in  the  force. 
The  instability  is  expected  to  occur  for  the  classical  motion  of  a  wide  variety  of 
molecules  with  more  than  one  oscillatory  mode. 

Our  results  are  important  for  two  reasons.  First,  the  classical  linear  motion  of  a 
rectilinear  symmetric  molecule  is  commonly  taught  in  physics  as  an  interesting  example 
of  linear  normal  mode  theory.  The  energy  of  each  normal  mode  remains  constant  in  this 
case.  We  have  shown  that  nonlinearity  can  dramatically  alter  this.  Second,  the  existence 
of  the  instability  raises  interesting  questions  about  the  oscillatory  behavior  of  actual 
(quantum  mechanical)  molecules.  As  example,  the  specific  heat  of  a  gas  of  the  molecules 
may  be  lowered.  As  another  example,  the  transition  probability  to  higher  energy  levels 
may  be  enhanced,  which  could  be  exploited  in  the  operation  of  some  lasers. 

The  instability  occurs  for  amplitudes  above  a  threshold  value,  where  this  value 
decreases  as  the  difference  in  frequency  of  the  two  normal  modes  is  reduced.  For  a 
softening  nonlinearity,  the  lower-frequency  normal  mode  can  be  unstable,  while  for  a 
hardening  nonlinearity,  the  upper- frequency  normal  mode  can  be  unstable.  Numerical 
simulations  agree  well  with  the  analytical  predictions,  although  some  interesting 
deviations  were  observed.  For  example,  parametric  resonances  other  than  the  principal 
resonance  can  occur  when  the  ratio  of  the  outer  masses  to  the  central  mass  is  not  small. 
In  addition,  if  the  mass  ratio  is  very  small,  the  principal  resonance  was  observed  to  occur 
when  the  theory  predicts  it  does  not  occur.  This  is  interesting  because  it  represents 
behavior  that  is  fundamentally  different  than  that  which  the  standard  leading-order 
perturbation  theory  predicts. 

I  constructed  an  apparatus  that  demonstrates  the  instability  in  the  case  of  both  the 
hardening  upper-frequency  mode  and  softening  lower-frequency  mode.  The  hardening 
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system  consists  of  gliders  on  an  air  track,  where  springs  and  repulsive  magnets  supply  the 
coupling.  In  accord  with  the  theory,  the  upper-frequency  mode  is  observed  to  be  unstable 
for  amplitudes  above  a  threshold  value.  The  resultant  motion  is  a  beating  between  the 
two  normal  modes.  Also  in  accord  with  theory,  the  lower-frequency  mode  is  observed  to 
be  stable  at  all  amplitudes.  Difficulty  in  observing  the  nonlinearity  in  the  hardening 
system  may  be  attributed  to  the  strong  linear  restoring  force  in  the  double  spring  coupling 
utilized  in  the  first  demonstration.  Additionally,  there  is  some  linearity  associated  with 
the  magnets.  Only  when  the  linearity  was  reduced  to  a  minimum,  the  nonlinearity  was 
increased  to  a  maximum  and  the  large  mass  was  increased  to  a  maximum  for  this  system 
was  the  instability  evident.  The  aforementioned  conditions  were  met  by  reducing  the 
system  to  a  single  spring,  extending  the  magnets  as  far  as  the  apparatus  would  allow,  and 
weighing  down  the  apparatus  as  much  as  the  system  could  handle  without  introducing 
friction. 

The  softening  system  consists  of  gliders  on  an  air  track,  where  slanted-springs  and 
attractive  magnets  supply  the  coupling.  Both  the  slanted-springs  and  the  attractive 
magnets  provide  a  softening  restoring  force  to  the  system.  In  accord  with  the  theory,  the 
lower-frequency  mode  is  observed  to  be  unstable  for  amplitudes  above  a  threshold  value. 
The  resultant  motion  is  a  beating  between  the  two  normal  modes.  Also  in  accord  with 
theory,  the  upper-frequency  mode  is  observed  to  be  stable  at  all  amplitudes.  Difficulty  in 
seeing  the  nonlinearity  in  the  softening  system  may  be  attributed  to  the  apparent 
weakness  of  the  nonlinearity  in  the  arrangement  and  a  quick  decay  of  the  initial 
amplitude  of  the  motion  supplied  to  the  system.  Even  after  the  first  slanted-spring  system 
was  optimized  to  supply  the  largest  softening  force  possible,  the  system  never  oscillated 
more  than  twice  above  the  amplitude  threshold  for  the  system.  Introducing  the  attractive 
magnets  into  the  system  helped  extend  the  lifetime  of  oscillations  above  the  threshold 
amplitude  so  that  the  system  could  exhibit  the  instability. 
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VI.  FUTURE  WORK 


According  to  our  extensive  literature  search,  an  instability  of  a  mode  due  to  the 
mode  parametrieally  driving  another  mode  has  not  been  investigated  for  quantum 
meehanieal  systems.  A  possible  interesting  eonsequenee  of  this  is  that  a  gas  of  triatomic 
moleeules  that  olassieally  exhibit  such  an  instability  may  have  a  speeifie  heat  that  is  less 
than  the  predieted  value.  This  should  be  analytically  investigated  and  experiments  eould 
be  performed. 

Another  interesting  quantum  meehanieal  possibility  that  should  be  investigated  is 
that  the  transition  probability  due  to  absorption  of  radiation  or  eollisions  eould  be 
enhaneed  in  eases  where  the  instability  cannot  proceed  beeause  it  would  violate  energy 
conservation.  This  could  have  a  benefieial  consequence  for  some  lasers. 

Our  numerieal  simulations  for  a  softening  nonlinearity  were  subjeet  to  the  “over- 
the-top”  instability,  whieh  oceurred  before  the  parametrie  instability  may  have  eeased  to 
oecur.  This  possible  upper  threshold  is  interesting  beeause  it  would  be  beyond  the 
standard  perturbation  theory.  The  effeet  ean  be  numerieally  probed  by  using  a  softening 
potential  energy  that  does  not  negatively  diverge.  An  example  of  sueh  a  potential  energy 
is  ln(l  +  X  )  (in  dimensionless  units). 

Our  numerieal  simulations  revealed  that  parametrie  resonanees  other  than  the 
prineipal  resonanee  readily  oeeur.  These  eases  ean  be  analytieally  investigated  and  the 
results  ean  be  compared  to  the  numerical  simulation  data. 

Our  successful  demonstrations  were  only  partially  quantitative;  we  did  not  fully 
quantify  them  in  order  to  compare  the  results  to  the  theory.  The  demonstrations  eould 
become  quantitative  experiments  if  the  force  vs.  displaeement  is  statieally  measured, 
whieh  would  allow  the  nonlinear  coefficient  to  be  determined.  Alternatively,  the 
nonlinear  eoefficient  could  be  determined  from  amplitude-dependent  period 
measurements.  The  nonlinear  coefficient  due  to  the  magnets  could  be  quantified  to 
ensure  that  the  demonstrations  exhibited  the  instability  due  to  the  principal  parametric 
resonance  and  not  due  to  nonuniformity  of  the  magnets. 
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I  did  not  create  a  demonstration  of  the  instability  for  a  triatomie  moleeule  with 
softening  nonlinearity,  using  only  slanted-springs  with  gliders  on  an  air  traek.  This 
arrangement  appears  to  be  inherently  only  weakly  nonlinear,  so  that  the  threshold 
amplitude  for  the  instability  is  diffieult  to  obtain  for  an  aetual  system.  More  work  needs 
to  be  done  with  the  slanted-spring  arrangement  to  aseertain  whether  it  ean  be  used  to 
convineingly  demonstrate  the  instability.  A  eomplieating  faetor  is  that  the  system  must 
be  symmetrie  so  that  simple  sets  of  initial  conditions  correspond  to  the  normal  modes  of 
the  system.  Otherwise,  linear  beating  will  occur,  which  would  mask  the  nonlinear 
beating  due  to  the  instability. 

Finally,  actual  molecules  do  not  have  symmetric  potentials.  That  is,  they  are 
approximately  modeled  by  both  eubic  and  quadratic  nonlinear  terms  in  the  foree.  The 
theory  for  this  should  be  done,  and  should  be  a  straightforward  generalization  of  our 
theory  with  just  a  eubic  nonlinearity.  The  theory  may  explain  why  our  preliminary 
asymmetrie  potential  arrangement  with  air  track  gliders,  springs,  and  magnets  failed  to 
exhibit  the  instability.  It  should  be  noted  that  an  alternative  demonstration  with  eubie  and 
quadratie  nonlinearities  ean  be  obtained  if  the  outer  or  inner  pairs  of  magnets  are 
removed  in  our  sueeessful  cubieally  nonlinear  glider  system.  It  should  also  be  noted  that 
this  demonstration  eould  beeome  a  quantitative  experiment. 
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APPENDIX  A:  NUMERICAL  SIMULATION  PROGRAMS 


A,  LOWER-FREQUENCY  MODE  PROGRAM 

lower.cpp 

This  program  determines  whether  the  lower-frequeney  (antisymmetrie)  mode  of  a 
triatomie  moleeule  is  stable  or  unstable  for  a  speeified  amplitude.  The  equations 
of  motion  are 

x"  +  X  =  alpha*(x^3  +  3*gamma^2*x*y^2) 

y"  +  gamma*y  =  alpha* (gamma^3*y^3  +  3*gamma*x^2*y)  , 

where  x  is  the  lower-frequeney  (antisymmetrie)  mode  and  y  is  the  upper- 
frequeney  (symmetrie)  mode. 


#inelude  <stdio.h> 
#inelude  <stdlib.h> 
#inelude  <math.h> 


#defme 

dt 

0.0001 

#defme 

epsilon 

0.0001 

#defme 

eriterion 

0.001 

#defme 

tmax 

100 

#defme 

alpha 

-1.0 

void  main() 

{ 

double  A,  rho,  gamma,  gamma2,  gamma3,  over; 

double  t,  X,  y,  u,  v; 

double  unonlin,  vnonlin; 

double  xsq,  ysq,  xnew,  ynew,  unew,  vnew; 

int  answer,  flag; 

MASSRATIO; 

printf("Enter  mass  ratio  m/M;  "); 
soanf("%lf',  &rho); 
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printf("\n"); 

gamma  =  1.0  +  2.0*rho; 
gamma2  =  gamma*  gamma; 
gammaS  =  gamma*gamma2; 

over  =  1.0  +  criterion; 

AMPLITUDE: 

printf("Enter  amplitude  of  lower  mode:  "); 

scanf("%lf',  &A); 

printf("\n"); 

ElEE  *fout; 

fout  =  fopen("dataout.dat",  "w"); 

X  =  A; 
u  =  0.0; 

y  =  epsilon; 
v  =  0.0; 

flag  =  0; 

for  (t  =  0.0;  t  <=  tmax;  t  =  t  +  dt) 

{ 

xsq  =  x*x; 
ysq  =  y*y; 

unonlin  =  alpha*(x*xsq  +  3.0*gamma2*x*ysq); 
vnonlin  =  alpha*  (gamma3*y*ysq  +  3.0*gamma*xsq*y); 

unew  =  u  +  (-X  +  unonlin)  *dt; 
vnew  =  V  +  (-gamma*y  +  vnonlin)*dt; 
xnew  =  X  +  unew*dt; 
ynew  =  y  +  vnew*dt; 

fprintf(fout,"%f  \t  %f  \t  %f  \n",  t,  x,  y); 

if(y  >  criterion)  flag  =  1; 

if  (alpha  >  0.0  &&  fabs(x)  >  over)  goto  OVERTOP; 

X  =  xnew; 
y  =  ynew; 
u  =  unew; 

V  =  vnew; 
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} 

fclose(fout); 

if  (flag  ==  1)  goto  UNSTABLE; 

QUESTION: 

printf("STABEE.  New  amplitude  (1),  new  mass  ratio  (2),  or  end  (0)?  "); 

seanf("%i",  &answer); 

printf("\n"); 

if  (answer  ==  0)  goto  END; 
if  (answer  ==  1)  goto  AMPEITUDE; 
if  (answer  ==  2)  goto  MASSRATIO; 
goto  QUESTION; 

UNSTABEE: 

printf("UNSTABLE:  Growth  of  upper  mode.  \n"); 
goto  QUESTION2; 

OVERTOP: 

printf("UNSTABEE:  Over-the-top  motion.  \n"); 

QUESTION2: 

printf("New  amplitude  (1),  new  mass  ratio  (2),  or  end  (0)?  "); 

seanf("%i",  &answer); 

printf("\n"); 

if  (answer  ==  0)  goto  END; 
if  (answer  ==  1)  goto  AMPEITUDE; 
if  (answer  ==  2)  goto  MASSRATIO; 
goto  QUESTION2; 


END: 


printf("End  of  program.\n\n"); 

} 
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B,  UPPER-FREQUENCY  MODE  PROGRAM 

upper.cpp 

This  program  determines  whether  the  upper-frequeney  (symmetrie) 
mode  of  a  triatomie  moleeule  is  stable  or  unstable  for  a  speeified 
amplitude.  The  equations  of  motion  are 

x"  +  X  =  alpha*(x^3  +  3*gamma^2*x*y^2) 

y"  +  gamma*y  =  alpha* (gamma^3*y^3  +  3*gamma*x^2*y)  , 

where  x  is  the  lower-frequeney  (antisymmetrie)  mode  and  y  is  the 
upper-frequeney  (symmetrie)  mode. 


#inelude  <stdio.h> 
#inelude  <stdlib.h> 
#inelude  <math.h> 


#defme 

dt 

0.001 

#defme 

epsilon 

0.0001 

#defme 

eriterion 

0.01 

#defme 

tmax 

1000 

#defme 

alpha 

-1.0 

void  main() 

{ 

double  A,  rho,  gamma,  gamma2,  gamma3,  over; 

double  t,  X,  y,  u,  v; 

double  unonlin,  vnonlin; 

double  xsq,  ysq,  xnew,  ynew,  unew,  vnew; 

int  answer,  flag; 

MASSRATIO; 

printf("Enter  mass  ratio  m/M;  "); 
soanf("%lf',  &rho); 
printf("\n"); 
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gamma  =  1.0  +  2.0*rho; 
gamma2  =  gamma*gamma; 
gammaS  =  gamma*gamma2; 

over  =  1.0/gamma  +  criterion; 

AMPLITUDE: 


printf("Enter  amplitude  of  upper  mode:  "); 

scanf("%lf',  &A); 

printf("\n"); 

ElEE  *fout; 

fout  =  fopen("dataout.dat",  "w"); 

X  =  epsilon; 
u  =  0.0; 

y  =  A; 
v  =  0.0; 

flag  =  0; 

for  (t  =  0.0;  t  <=  tmax;  t  =  t  +  dt) 

{ 

xsq  =  x*x; 
ysq  =  y*y; 

unonlin  =  alpha*(x*xsq  +  3.0*gamma2*x*ysq); 
vnonlin  =  alpha*  (gamma3*y*ysq  +  3.0*gamma*xsq*y); 


/* 

*/ 


unew  =  u  +  (-X  +  unonlin)  *dt; 
vnew  =  V  +  (-gamma*y  +  vnonlin)*dt; 
xnew  =  X  +  unew*dt; 
ynew  =  y  +  vnew*dt; 

fprintf(fout,"%f  \t  %f  \t  %f  \n",  t,  x,  y); 

if  (x  >  criterion)  flag  =  1; 

if  (alpha  >  0.0  &&  fabs(y)  >  over)  goto  OVERTOP; 


} 


X  =  xnew; 
y  =  ynew; 
u  =  unew; 
V  =  vnew; 
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fclose(fout); 

if  (flag  ==  1)  goto  UNSTABLE; 

QUESTION: 

printf("STABEE.  New  amplitude  (1),  new  mass  ratio  (2),  or  end  (0)?  "); 

seanf("%i",  &answer); 

printf("\n"); 

if  (answer  ==  0)  goto  END; 
if  (answer  ==  1)  goto  AMPEITUDE; 
if  (answer  ==  2)  goto  MASSRATIO; 
goto  QUESTION; 

UNSTABEE: 

printf("UNSTABEE:  Growth  of  lower  mode.  \n"); 
goto  QUESTION2; 

OVERTOP: 

printf("UNSTABEE:  Over-the-top  motion.  \n"); 

QUESTION2: 

printf("New  amplitude  (1),  new  mass  ratio  (2),  or  end  (0)?  "); 

seanf("%i",  &answer); 

printf("\n"); 

if  (answer  ==  0)  goto  END; 
if  (answer  ==  1)  goto  AMPEITUDE; 
if  (answer  ==  2)  goto  MASSRATIO; 
goto  QUESTION2; 


END: 


printf("End  of  program.\n\n"); 

} 
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APPENDIX  B:  HIGHER  ORDER  PARAMETRIC  RESONANCES  IN 
LOWER  FREQUENCY  MODE  WITH  HARDENING 
NONLINEARITY 


The  following  graphs  demonstrate  the  different  parametrie  resonanees  and  ean  be 
understood  by  examining  the  ratio,  h/fd,  of  the  response  frequeney  to  the  drive  frequeney. 
For  the  prineipal,  isoehronous,  third  order,  and  fourth  order  parametrie  resonanees,  h/fd  = 
14,  1,  3/2,  2,  respeetively.  Due  to  the  nonlinearity,  the  mode  being  parametrieally  exeited 
experienees  a  drive  at  twiee  the  frequeney  of  the  initial  mode. 


Figure  B.  1 .  Example  of  principal  parametric  resonance  in  the  lower- frequency  hardening  mode,  p  = 
0.009,  A  =5.26300,  6t  =  0.001,  epsilon  =  0.0001,  criterion  0.001. 


Figure  B.2.  Example  of  isochronous  parametric  resonance  in  the  lower-frequency  hardening  mode,  p 
=  0.5,  A  =8.2693,  5t  =  0.001,  epsilon  =  0.0001,  criterion  0.001. 
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Figure  B.3. 


Figure  B.4. 


Example  of  third  order  parametric  resonance  in  the  lower-frequency  hardening  mode,  r  = 
2.0,  A  =  94.09,  5t  =  0.0001,  epsilon  =  0.0001,  criterion  0.001. 


Example  of  fourth  order  parametric  resonance  in  the  lower-frequency  hardening  mode,  r 
=  4.2,  A  =  6.006,  5t  =  0.0001,  epsilon  =  0.0001,  criterion  0.001. 
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APPENDIX  C:  NONLINEAR  BEHAVIOR  OF  A  SLANTED-SPRING 

OSCILLATOR 


A  demonstration  of  a  classical  triatomic  molecule  in  which  the  coupling  hardens 
(slope  of  force  vs.  displacement  increases  with  amplitude)  can  readily  be  obtained  with 
springs  and  magnets  (Sec.  IV.D).  We  also  desire  a  demonstration  in  which  the  coupling 
softens  (slope  of  force  vs.  displacement  curve  decreases).  Although  it  is  not  obvious, 
such  a  situation  can  be  obtained  with  a  slanted-spring  arrangement  (Fig.  C.l).  In  this 
section,  we  determine  the  leading  (cubic)  order  nonlinearity  in  terms  of  the  parameters  of 
the  system,  and  thereby  determine  the  conditions  for  which  the  slanted  spring  oscillator 
softens. 


Figure  C.l.  Slanted-spring  oscillator  in  equilibrium.  The  mass  oscillates  back  and  forth  as 
denoted  by  the  arrow.  This  is  a  nonlinear  oscillator;  Flooke’s  law  only  holds  approximately  for  small 

amplitudes. 

We  let  the  unstretched  length  of  each  spring  be  c.  For  any  displacement  x,  we  can 
determine  the  force  from  Fig.  C.2.  For  practical  purposes,  it  is  assumed  that  |x|  <  a  and  c 
is  sufficiently  small  such  that  each  spring  never  compresses.  However,  the  theory  holds 
for  any  displacement  and  for  any  values  of  the  parameters.  The  net  force  on  the  mass  m 
is 


F  =  -2k 
+  2k 


(a  -I-  x)  -b  -  c  cos(0) 


a  -  x)  -I-  -  c  cos((p) 


(C.l) 
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Substituting  the  relationships 


/n\  a  +  x  .  /  X  a-x 

cos(9)=  I  ^  =  and  cos((p)  = 


(a  +  x)  +b 


(a-xf +b^ 


and  simplifying,  we  find  that 


F  =  -2k 


+  2k 


1- 


1- 


(a  +  x)  +b' 


(a-x)  +b^ 


(a  +  x) 


(a-x) 


(C.2) 


This  gives  the  exact  force  as  a  function  of  the  displacement. 


Figure  C.2.  Essential  geometry  of  a  slanted-spring  oseillator  displaced  a  distance  x  from 

equilibrium. 


Before  expanding  Eq.  (C.2)  in  a  Taylor  expansion  in  x,  we  consider  a  graph  of  |F| 
vs.  X  for  k  =  1/4  and  a  =  b  =  c  =  1 .  The  result  is  shown  in  Fig.  C.3,  which  reveals  that  the 
force  softens.  This  is  surprising.  As  the  displacement  increases,  the  restoring  springs 
become  more  longitudinal  and  the  anti-restoring  springs  to  become  more  perpendicular, 
which  causes  the  slope  of  force  vs.  displacement  to  increase.  That  is,  in  Eq.  (C.l)  the 
cos(0)  factor  increases  and  the  cos((p)  factor  decreases.  Hence,  it  is  expected  that  the 
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system  hardens.  However,  in  addition  to  this  effeet,  the  length  of  each  spring  softens 
with  increasing  x.  This  second  effect  evidently  dominates  the  first. 


Fig.  C.3.  Exact  (nonperturbative)  magnitude  of  the  force  vs.  displacement  for  a  slanted-spring 
oscillator.  The  parameters  are  k  =  1/4  and  a  =  b  =  c  =  1,  where  c  is  the  unstretched  length  of  each  spring. 
The  dashed  line,  which  has  slope  1  -  is  Flooke’s  law  in  the  limit  of  small  displacements.  The  dotted 
curve  is  the  result  of  the  perturbative  theory,  which  is  valid  to  third  order  in  the  displacement. 


For  small  displacements,  we  can  expand  Eq.  (C.2)  to  third  order  in  x.  This 
requires  a  substantial  amount  of  algebra.  We  begin  by  noting  that 


(a  -b  x)^  -b  -bb^.  1  -b 


2ax  -b  X" 


a"  +  b" 


(C.3) 


2  2  1/2 

The  quantity  [(a  -  x)  +  b  ]  can  be  obtained  by  simply  replacing  x  with  -x  in 
this  expression.  Next,  we  use  the  expansion 

1  ,1  3263 

I -  =  1 - £-b-£^ - £'^-b---  . 

Vi+i  2  8  16 


2  2  2 

Substituting  £  =  (2ax  +  x  )/(a  +  b  ),  retaining  terms  only  to  cubic  order  in  x, 
grouping  terms  according  to  the  power  of  x,  and  simplifying,  yields  the  following  from 
Eq.  (C.3): 
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Va"  +  b"  ,  a  2a" -b"  2  2a' -Sab" 

P  +  xf  +  b^  a^  +  b^  2(a^  +  b^f  2(a^  +  b^f 


In  view  of  Eq.  (C.2),  we  multiply  this  by  a  +  x,  expand  to  third  order  in  x,  and 
simplify.  The  result  is 


Va^ 
^[a  +  xf  +b" 


(a  +  x) 


=  a  + 


b"  Sab"  2  4a"b"-b" 

^5 - 5-X - ^X  + - ^X 

a  +b^  2(a"+b")  2(a"  +  b") 


Finally,  we  substitute  this  expression  into  Eq.  (C.2)  to  determine  the  force.  Due 
to  the  symmetry,  all  even  terms  in  x  cancel,  and  the  result  is 


F  =  -4k 


1- 


b"c 


(a"  +  b"') 


3/2 


x  +  2k 


4a" -b" 

(a"  +  b"p 


b"cx' 


(C.4) 


There  are  several  special  cases  of  the  result  (C.4)  that  can  be  considered.  First, 
note  that  if  b  =  0  (longitudinal  oscillator),  then  the  equation  yields  the  correct  result  F  =  - 
4kx,  independent  of  c.  Second,  if  a  =  0  (transverse  oscillator),  then  Eq.  (C.4)  reveals  that 
the  oscillator  is  always  hardening  if  c  0.  Third,  if  c  =  0,  then  Eq.  (C.4)  yields  F  =  -4kx, 
independent  of  a  and  b.  This  is  perhaps  a  surprising  result,  although  from  Eq.  (C.2)  it  is 
immediately  seen  to  be  valid  exactly  (to  all  orders  of  x). 

As  shown  in  Fig.  C.3,  for  the  parameters  k  =  1/4  and  a  =  b  =  c  =  1,  the  result  (C.4) 
agrees  very  well  with  the  exact  force  if  the  displacement  is  small.  The  deviation  is  a 
maximum  at  the  maximum  displacement  of  unity,  as  expected,  and  this  deviation  is  7.5%. 
The  result  (C.4)  shows  that  the  slanted-spring  oscillator  softens  if  2a  >  b,  provided  b  0 
and  c  0.  We  confirmed  this  approximate  result  by  plotting  the  exact  force  (C.2)  for  k  = 
1/4  and  a  =  c  =  1,  with  b  =  1.8  and  2.2.  The  curves  are  nearly  straight  lines,  with  the  b  = 
1.8  curve  having  negative  curvature  and  the  b  =  2.2  curve  having  positive  curvature. 
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